Pulsar scintillation patterns and strangelets 
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We propose that intersteUar extreme scattering events, usuaUy observed as pulsar scintillations, 
may be caused by a coherent agent rather than the usually assumed turbulence of H2 clouds. We 
find that the penetration of a flux of ionizing, positively charged strangelets or quark nuggets into a 
dense interstellar hydrogen cloud may produce ionization trails. Enhanced electron densities would 
then form and constitute a lens-like scattering screen for radio pulsars and quasars. 



A variety of scintillation phenomena observed from 
pulsars and quasars require interstellar scattering screens 
that contain compact regions of high electron density. 
These include quasar Extreme Scattering Events (ESE) 
UIj 01 , pulsar parabolic arcs [j, |4 1 and Galactic Center 
scattering of OH maser sources 0, Q. Many of these 
phenomena, in particular the ESEs, require enhanced 
electron density regions of A.U. (~ 1.510"'^'^ cm) size. 
The overpressure P/ks ~ 10^ — 10* K cm~^, as estimated 
from typical temperatures T ~ 10^ K and particle densi- 
ties n ^ 10^ — 10'' cm~'^, is difficult to explain in any con- 
ventional scattering screens embedded in dense molecular 
H2 clouds. The only plausible environment where such 
pressures might be attained would be in the dense cloud 
cores (cf. the compact ionized cloud model developed by 
Walker [71, |8[ ) . The required source properties in this lat- 
ter model require a significant mass in extremely dense 
cold gas clumps to source the ionized gas. The stabil- 
ity of such cold clumps is questionable, although exotic 
models have been proposed [9|. 

There is some direct evidence that the ionized clouds 
arc highly elongated [10[. It is an already well-known 
fact that plasma lenses result from electron (over-) under- 
densities. Electron over-densities result in a faster phase 
velocity, corresponding to a concave (divergent) optical 
lens, while under-densities are associated with a conver- 
gent lens ll[. Quantitative lens models for the ESEs [7[ 
estimate column densities pe ~ 10^^ cm^^. 

Here we are interested in estimating the effect of the 
enhancement of the electron column density on A.U. 
scales as a result of the formation of ionized trails in HI 
and H2 clouds. We propose a mechanism that provides an 
alternative to postulating the existence of controversial 
clumps of dense molecular hydrogen, and naturally gen- 
erates pervasive ionized trails in dense interstellar cloud 
cores. Our model invokes strangelets |12| (also known as 
nuclearites) , finite droplets of quark matter with non-zero 
strangeness fraction and slightly charged. They are cur- 
rently being searched on earth [13[ for as final products 
in heavy ion collisions with the ALICE experiment at the 
LHC, in the CDMSII under the form of light ionizing par- 
ticles or in the AMS-02 mission. It is expected that these 



quark droplets can be naturally generated by a series of 
different astrophysical events where a nucleon-quark de- 
confinement transition may take place, e.g. neutron star 
(NS) collisions, NS or black hole combined binary merg- 
ers or in NS to quark star (QS) conversions. This latter 
process may be induced by internal heating due to dark 
matter (DM) annihilations [1J| (under the assumption of 
a Majorana DM particle candidate) and even leave ob- 
servable traces in the pulsar distribution [l5| or in the 
emission of very short gamma ray bursts (GRBs) with 
typical time scales T90 ^ 0.1 s [16|- Due to the large 
gravitational and nuclear binding energies released in the 
transitioning process, a mass ejection episode is expected, 
possibly seeding the interstellar medium with a fraction 
of strangeness-carrying lumps of matter formed during 
the phase of nucleon deconfinement. Energetics for short 
GRBs isotropic equivalent photon emission show that the 
measured value E^igo ~ 10^* — 10^^ erg is compatible 
with relativistic mass ejecta A/oj ^ 10"'' M© able to con- 
sistently produce observable gamma rays. It has been 
actually proposed that quark matter droplets might par- 
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tially populate cosmic ray (CR) primaries (see e.g 

In this context, let us consider a volume of interstellar 
space that initially contains a '-^ 30 pc cloud radius of H2 
(molecules) and HI (neutral) atoms. Typical ionization 
reactions are of the form, X -I- H2 — ?> X -t- H2' + e" , or 
X -t- HI ^> X + H+ + e~ , where X is the incoming charged 
strangelet (assumed to have unit charge in the previous 
reaction, but more generally with charge Ze) . The energy 
needed to ionize a hydrogen atom (molecule) , initially in 
the ground state is I(HI) = 13.6 eV (I(H2) = 15.6 eV). 
In astrophysical CGS units we use the conversion fac- 
tor leV = 1.62 X 10~^^ erg. The dimensions of the 
cloud can be estimated assuming a radius for the cloud 
core Re ~ lOpc ^ 3 x lO^'^ cm, then the core vol- 
ume is Vc w 47rRy3 ~ 1.1 x 10^^ cm^. Since the 
atomic number density in the cloud core is uq ^ 10^ 
cm~'^, the number of particle species in the original neu- 
tral cloud (HI, protons, electrons) is accordingly given 
by TiQ Vc — 1.1 X 10^"^. In general, other ionizing agents, 
such as CRs can cause ionization of H2 or HI with an 
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FIG. 1. Physical scenario of ionization in the H2 cloud due 
to an emitting source of strangelets in the line of sight to the 
observer. 
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A more exotic flux of potentially ionizing electri- 
cally charged strangelets through the cloud will depend, 
on the one hand, on the time and spatial distribution of 
their astrophysical sources (besides a possible primordial 
background) and, on the other hand, on their peculiar 
nature. 

An estimate of the incoming flux starts with the 
strangelet formation rate associated with the different 
astrophysical sources involved. We will consider a sim- 
plified model with a unique A baryonic number lump 



of mass 



'STR 



< Attin, where ttt-n is the nucleon mass. 



The strangelet number production rate in the astrophys- 
ical ith-process will be given by dN\/dt = 77iAfoj,i/m.sTR- 
Mej^i is the mass rate ejected in the ith-process and rji is 
the efficiency of strangelet ejection for the ith-process in- 
volved. For example, for NS collisions it is expected that 



Mej '- (10-5 - 10-^)Mq while in a NS merger the mass 
ejection is Moj '^ (5 x 10^^ — 7 x lO-'^)Af0 for equal-mass 
binaries with total mass m — 2.7Mq [2^. In a NS to QS 
transition capable of emitting a GRB, it is expected that 
Afcj < IO-'^A'/q. According to [l6| the rate of these tran- 
sitions is i?NS-QS — (8 X IQ-^ — 3 X 10-^)i?sNtii being 
i?SNtii — 10^^ yr^i the rate of type II supernovae in our 
galaxy. 

Having enumerated the possible processes that may 
constitute a source of the strangelet flux, for the rest 
of this work wc will consider a generic source where 
the dcconflncmcnt transition can take place with galac- 
tic appearance rate R = R c:^ IQ-^ yr-^ and ejected 



mass Mf. 



M, 



10 A/0. The ejected mass rate 
is then Afej = RM^j ~ lO-^^Afoyr-i. Since the effi- 
ciency of strangelet production depends on the so far un- 
known details of the engine model, we will consider that 
only a small fraction 77 ^ 10-^ is ejected under exotic 
form [16|. Nevertheless, strangelets should be emitted 
with A- values larger than a critical stability value J23l |. 
A > ^min — 10"'^ — 10^ so that they can possibly decay to 
the lightest energetically stable ^min-species. 



The strangelet number production rate at the source 



dNA 
dt 



2x10"^ 




Mr 



10-10 A/0 yr-1 



(1) 



It is expected that a possible emission distribution func- 
tion at the source fs{Z,f3) can modulate this rate. We 
will not consider this refinement here, and in what fol- 
lows we will assume fs{Z^P) = 1. From such a source, 
the corresponding unscreened fiux measured by an ob- 
server at distance dso ^ 600 pc is given by 



F 



dNA 



1 



dt 47r(i|Q 



1.6 X 10"" A 



-1 _9 _1 _1 

cm sr s , 



(2) 



^l3 



for A > ^min- Note that for strangelets being currently 
searched for in neutrino telescopes on earth, there is a 
lower limit A > lO^"^ GeV/c^. Our estimates yield in 
that case F < 1.6 x lO^^^ cm^^sr-is"!. Unfortunately, 
this is not presently within the reach of ANTARES, who 
report[2J| a testing capability of 90% C. L. flux limit of 
-Fantares ^ 2 X 10" 
lccCubc-22, who report F1C22 
A> 10^^ nuclearites ^. 

Once the nuclearites are produced it may happen that 
their trajectories intersect with a II2 cloud along the line 
of sight to the observer, as depicted in Fig. [1] In this 
scenario one can estimate A^^ , the number of ionizing 
A— sized strangelets traversing the cloud in its lifetime, 

T, ^ 10^ yr, as Nf = m^^ Jo' ^dt. Using typical 
values [26| for the distances in this expression, we obtain 
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Nf ~ 1.25 10'*' 



A- 



0.01 



Re 
Ipc 



200 pc 

dsc 



Wyr 



(3) 
where dsc is the distance from the source to the cloud. 
Since beaming is observed in some very short GRBs (ex- 
pected in this scenario |l6j) with opening angle 9-^ ^ 

1° — 30°, an efficiency factor yyn ~ -J- has been intro- 
duced. 

The peculiar nature of the quark droplets is highly un- 
certain but it is usually assumed that their charge is small 
and distributed positively on the surface [l7[. For ordi- 
nary strangelets Z = 0.3^^/^ [l^ while for CFL (color- 
flavour- locked) strangelets Z ~ Q.iA^I^ . Even smaller 
charge-to-mass ratios are energetically possible for in- 
termediate masses. A ^ 10^ — lO^®, assuming a strong 
coupling as ~ 0.9 |28|. For example, for A ~ 10^, 
Z/A ~ lO"'' while larger strangelets A ~ 10^^ have 
Z/A ~ 10"^ - 10-3 and even Z/A < 0. In this work 
we constrain strangelets to have Z > 1. 

As the slightly charged and heavy strangelet traverses 
the hydrogen cloud in nearly straight lines, it continu- 
ously loses kinetic energy. Typically, the ejection energy 
at the source allows these particles to achieve a Lorcntz 
factor bounded by a saturation value, 7 ^ 7sat with 
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FIG. 2. Logarithm of the kinetic energy loss rate in the cloud 
as a function of the logarithm of the incoming {A = 10'') 
strangelet kinetic energy per baryon number. Several charge 
states are shown and a CR proton case with Z = 1 is depicted 
for comparison. 
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FIG. 3. Logarithm of the strangelet range in the cloud as a 
function of the logarithm of the incoming kinetic energy per 
baryon number. We consider a CR proton Z = 1, ordinary 
and CFL strangelets with A — 10^ and A — 10^ respectively, 
and larger strangelets with A — 10^^ and A — 10^* with 
Z/A = 10"^ 



7sat '^ 20 — 1000 [16|. Strangelets of net effective charge 
Zi may alter the cloud ionization fraction either by ion- 
izing hydrogen or gaining electrons so that they change 
its own incident state of charge, as measured for massive 
ions on gas [29|. For low velocities (/? = v/c) this may be 
important as it diminishes the ionizing power from the 
bare charge Z to an effective Zi = Z{1 - e-(0-95/^/"^''") 

Interaction in the HI cloud may arise due to a variety 
of processes [28| jSlI . The strangelet (kinetic) energy loss 
per unit length due to ionization and pion production 
can be obtained from the stopping power ^ = ^4^- At 
lower energies, ionization is the most important process, 
and for an incoming particle with effective charge Zi and 
velocity /3 the energy loss rate reads [32| 

f = -1-82 X 10-3 (_^) [ZMZ, PHEiAP)] eV/s 



10 '*pc ^, assuming 



10« 



(4) 



where 



^iZ,(3) = 



[1 + 0.0185 ln(/?)6l(^-^o)] 



/33 + 2/33 



and 



5(A,/3) = 0.72A°-53^i-286i(7 



1.3), 



6{x) is the Heaviside function, 7"^ ~ ^\ — fP- and 
/3o = 0.01 is the electron orbital velocity. This func- 
tion is shown in Fig. [2] versus the strangelet kinetic en- 
ergy per baryon number for an A = 10^ case. Different 
amounts of charge ZjA = 5 lO^^ (sohd), ordinary (dot- 
ted) and CFL (dot-dashed line) strangelets and a CR 
proton case (dashed line) have been considered. For rel- 
ativistic strangelets most of the energy is kinetic and they 
show an enhanced ionizing power with respect to the lin- 
ear case and a change in slope around TjA ^ 500 MeV 
due to pion production for larger energies. 

To estimate the hydrogen ionization rate by 
strangelets, we consider the number density of old NSs 



in our galaxy to be ?ioidNS 

NSs in a sphere of radius ~ 10 kpc. The number den- 
sity of molecular clouds (MC) is n^c ^ lO^^pc^^^ 
so we require on average that a fraction /conv '^ 1% 
of old NSs convert to quark stars and be sources of 
ionizing strangelets. The strangelet injection rate can 
be obtained using the production rate per source from 

Eq. (HD, dNA/dtUoMNsfconv ^ 2 X W k^^^T-^^^-^ = 

2.3 X 10-24A-icm-3s-i. 

Taking an escape time from the MC by diffusion of 
the -few GeV/A strangelet as At - Z^I'^Wyx^ 10^° 
s, the strangelet number density is tt-str = 2.3 x 
10-i42-i/3^-icm-3. This is much lower than the CR 
density at the low-energy break in the spectrum of Galac- 
tic CRs at '^ 1 GeV, computed from the CR flux For ~ 
lO-^m-^s-isr-i orncR-310-i°cm-3 [s^. 

The hydrogen ionization rate C^ averaged over 

the molecular cloud, is now estimated to be (^ — 

{dT^°^^/dt) nsTR{I)~^nJi^, per H atom assuming (/) ~ 36 

2/3^, (5) eV per H ionization. Then C^ = ^^1.2 x lO'^^ g-i, per 

HI atom for a strangelet with A = 10^. 

The strangelet mean free path or average range can 

be calculated as i?sTR — I ^"^ i^) ■ ^^ strangelets 
have a net charge Z then the new range will be smaller 
since i?sTR = Rstr{z=i)/Z^, and therefore will provide 
a more efficient ionization. In Fig. [3] we show the loga- 
rithm of the strangelet range in the cloud as a function of 
the logarithm of the incoming kinetic energy per baryon 
number. We consider a proton CR case Z = 1, ordinary 
and CFL strangelets with A = 10^ and A = 10^ respec- 
tively, and larger strangelets with A ~ 10^^ and A = 10^^ 
with Z/A = 10"^. As seen the cloud can effectively stop 
the very heavy nuggets while this is not possible for the 
smaller droplets even at lower kinetic energies. 

Once the cloud core is populated with trails of Hll 
(ionized hydrogen), this region contains H"*" ions that 
tend to recombine radiatively as H+ -I- e~ — )■ H -I- hi/. 
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FIG. 4. Logarithm of the electron column density in the cloud 
as a function of the logarithm of the incoming kinetic energy 
per baryon number. We consider a CR proton Z = 1, ordinary 
strangelet with A = 10^, and large strangelets with A = 10^ 
and Z/A = 5 x 10"*, A = lO" and Z/A = 10"=* and A = 10^^ 
and Z/A = 10"*. 



The balance between ionization and recombination de- 
termines the ionized fraction. To estimate the possibility 
of a net ionization in the HI cloud one must compare 
the ionization and recombination times, Tjoniz and Tioc 



< 



If recombination is 



respectively, so that Tic 
slower, there would be a net amount of ionized HI in 
the cloud to cause an electron density enhancement. In 
addition, there might be a fraction of net charge due to 
the accumulation of positively charged lumps stopped by 
the cloud. The recombination time, t^oc, is calculated 
in the on the spot approximation as t^oc = l/necS^' = 
3.85 X lQ^'^n-^T2-^ s where T4 = T/10'' K. We estimate 
the average ionization time from the average number 
of ionization in the cloud lifetime as rioniz — Tc/N'^. 
Since iioniz << ^rcc then it is indeed possible to ob- 
tain net ionization. The ionization of interest for the 
ESEs occurs along the trajectory, that spans the MC, 
but only over a trajectory width given by the prod- 
uct of hydrogen recombination time and sound speed, 



10^ cms ^, namely h ~ c^a 



H "-H 



10 



15^ 



cm. 



Therefore the local ionization rate is increased with re- 
spect to a Z = 1 CR by ~ Z^{Rc/h)^ - lO^Z^nl 



to give C^ - Z^ lO-i^s" 



M' 



per HI atom. Equating 



4 X 10- 



with (h yields an electron density of Ue 



Z n\cTa~'^ . Typical events have been argued 
to be explained by electron densities when the MC gas 
will be fully ionized in a small volume around the track 
of about [3JI rif, = (100 — 400) cm"'^, with some cases up 
to rig = lO^'^cm"'^. Assuming an A = 10'^ nugget and i. 
e. Z = 1 — 30 depending of its nature, this yields trail 
widths of /i - 3.3 x(10-5-10-6)i?c -^3 x 10^3-3 x lO^^ 
cm as required by the ESE events. 

The electron column density, pe, can be calculated 
along the line of sight (los) as pe = /j^g ne{s)ds, then the 
enhancement coming from the ionization produced by the 
incoming particles is estimated as Ape ~ /j^^, Aueds. In 
Fig. |3]we show the logarithm of the electron column den- 
sity in the cloud as a function of the logarithm of the in- 



coming kinetic energy per baryon number. We consider a 
CR proton Z = 1 (dashed line) , ordinary strangelet with 
A = 10^ (dotted), and large strangelets with A = 10^ 
and Z/A = 5 x 10"^ (short dashed), A = 10^^ and 
Z/A = 10-3 (sohd) and A = lO^^ and Z/A = 10"^ (dot- 
dashed). We have used the ionization produced with A''^* 
particles, the more heavy and charged the more effective 
are the particles at causing scintillation. These values are 
compatible with some quantitative models of ESE [8| . 

We have shown that pulsar scintillations may be caused 
by an ionization agent constituted by positively charged 
lumps of strange matter. Typical obtained over-densities 
in electron column densities are compatible with those 
of quantitative ESE models. Compact A. U. sized re- 
gions are involved in the geometry of these events. Since 
these quark nuggets remain to be experimentally discov- 
ered further work is needed to experimentally confirm 
this theoretical scenario. 
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